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Abstract 

Using the experimental data from the STAR, PHENIX, ALICE and CMS programs on the 
rapidity and energy dependence of the pT spectra in p+p collisions, we show that a universal 
distribution exists. The energy dependence of temperature T and parameter n of the Tsallis 
distribution are also discussed in detail. A cascade particle production mechanism in p+p collisions 
is proposed. 
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I. INTRODUCTION 


The particle spectrum is a basic quantity measured in experiments and it can reveal the 
information of particle production mechanism in heavy-ion collisions. Recently, the Tsallis 
distribntion has attracted many theorists’ and experimentalists’ attention in high energy 
heavy-ion collisions jl-24|. The excellent ability to fit the spectra of identified hadrons and 


charged particles in a large range of pt up to 200 GeV is quite impressive 


20 


231. From 


the phenomenological view, there may be real physics behinc 


ogy work, e.g. Regge trajectory for particle classification 


the prominent phenomenol- 


25| . p-|-p collision experiments 


have been performed and measured under different energies. Since p-|-p collision is very 
simple compared to nucleus-nucleus collision, the measurements of p-|-p collisions are used 
to understand the particle interaction, particle production mechanism and as a baseline for 
nucleus-nncleus collisions. Many efforts have been pnt to stndy the particle spectra prodnced 


in p-|-p collisions using Tsai 
used in the literatnre [l 
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is distribution. Different versions of the Tsallis distribution are 
. The parameter T in the Tsallis distribntion was interpreted 


as temperature. All of them can fit the particle spectra very well but they give different 
temperatnres. In this work, we wonld like to stndy the connections and differences among 
different versions of the Tsallis distribntion. We collected p-|-p collisions data with different 
Pt ranges and different rapidity cuts from different experiment groups at RHIC and LHC 
and did a systematic study of the particle spectra using one of the Tsallis distribntions. 

The paper is organized as following. In section II, we review different versions of the 
Tsallis distribntion used in the literatnre and analyze their connections and differences. We 
also give the form of Tsallis distribution used in our anlysis. In section III, we show our 
results of particle spectra from p-|-p collisions and analyze them. We compare our results 
with the ones obtained by the other anthors in the literatnre. A brief conclusion is given in 
the section IV. 
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II. TSALLIS DISTRIBUTIONS 


The STAR fl, PHENIX |2[ Collaborations at RHIC and ALICE 
Collaborations at LHC adopted the form of Tsallis distribution 


5| and CMS 6| 


dp^ 


d‘^N 


2tipt dydpT 
dN [n — 1)(?7, — 2) 




( 1 ) 


dy 27rnC[nC + m{n — 2))] 

where rriT = is the transverse mass, m was used as a fitting parameter in ref. [l|, 

but it represents the rest mass of the particle studied in refs. [j-IgJ. n and C are fitting 
parameters. With this form of the Tsallis distribution, when pt ^ m, we can ignore the m 
in the last term in Eq. ([T]) and obtain oc This result is well known because high 

energy particles come from hard scattering and they follow a power law distribution with 


Pt- When px ^ m which is the non-relativistic limit, we obtain rriT ~ ^ = Et 

frclassical 

and oc e”— 

dp'^ 


classical 


c , i.e. a thermal distribution. The parameter C in Eq. ([T]) plays the 
same role as temperature T. 


In refs. 


14 


18| . the following Tsallis form is used 


d^N mrcoshp mrcoshp-/! ^ 

= gV [1 + {q- 1)-H-9. 


dp^ 


(2vr) 


T 


( 2 ) 


based on thermodynamic consistency arguments. Where g is the degeneracy of the particle, 
V is the volume, y is the rapidity, p is the chemical potential, T is the temperature and q 
is the entropic factor, which measures the nonadditivity of the entropy. In Eq. ([2]), there 


are four parameters V, p, T, q. p was assumed to be 0 in refs. 14j-ll7l| which is a reasonable 


assumption because the energy is high enough and the chemical potential is small compared 
to temperature. In the mid-rapidity y = 0 region, Eq. (j2]) is reduced to 


dp^ 




{271 


T 


( 3 ) 


It becomes very similar to Eq. CDi, but there are some differences, i.e. rriT replaces rriT — m 
in the bracket and there is a term rriT in front of the bracket as well. We have seen that two 
kinds of representation have been used: one is with parameter n, i.e. Eq. ([1]) and the other 
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is with parameter q, i.e. Eq. ([3]). There is no direct match between n and q because Eqs. 
([H |3]) are just similar. We find a connection between n and q in the limit at large pt- 
When Pt m, from Eq. (|3]), we can obtain 


dp^ 


oc Pj, 


1 

q-1 


(4) 


Recalling that oc when pt ^ th from Eq. ([T]), therefore the relation between n 

and q is 


n = 


q - 1 


(5) 


Another treatment to find the relation between n and q can be found in ref. 241. 

In ref. Sena et al. applied the non-extensive formalism to obtain the probability of 
particle with mometum pt 


1 da 
a dpT 


cpT j dpL[l + (g - 1)13 \Jp\ +Pt + 


( 6 ) 


where c is the normalization constant, g is a parameter, (^ = -f and m is the mass of particle. 
With the approximation px very large compared to pi and m 26(], Eq. ([ 6 ]) can be rewritten 
as 


1 da 
a dpT 


c[2(g-l)]-V2i?(i,^ 


+ (g- l)n] 



(7) 


where u = ^ and B(x, y) is the Beta-function. We can see that Eq. ((Tj) is also similar to 
Eqs. dDEl) but not exactly the same. We repeat the same limit condition, let px be very 
large, g 7 ^ 1 and (g — l)u ^ 1 , then 


1 da 
apxdpx 


_ ^ 

(xpx''~\ 


( 8 ) 


which is the same as Eq. (jl]). 

In the previous versions of the Tsallis distribution, the rapidity cuts for the experimental 


data are not taken into account. In ref. 


2l|, Wong et al. proposed a new form of the Tsallis 
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Figure 1: (Color online) Fitting results using the Tsallis distribution Eq. 
and p in p+p collisions at y/s = 62.4 GeV. Data are taken from PHENIX 


^^i^^ or 7r°, 


TT^, p 


distribution function to take into account the rapidity cut, 


dpS IliKo - 


, dy.d^N 


(9) 


Where 


^iV,PT) = 
dp 


' 1 - 




2 


cosh y 


( 10 ) 


with 


and 


y = - In 
^ 2 


■\J cosh^ p + m? + pt sinh p - 
\Jp\ cosh^ p + m? — Pt sinh p 


^ = r^n.^)-n 
dp^ dy^ riT’ ’ 


Erp = rriT — m, 


where is assumed to be a constant in ref. 

dy 

We can rewrite Eq. CD into 


2 ll. 


( 11 ) 


^ + EZ)-r 

dp^ dp dy nT' 


( 12 ) 


From the experimental measurements, we know that ^ is almost a constant when p is not 
large j^. We can simply treat it as a constant. Subsituting Eq. flT^ into Eq. (El), we 
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Figure 2: (Color online) Fitting results using the Tsallis distribution Eq. ()13p for p 

and p in p+p collisions at ^/s = 200 GeV. The data in (a) are taken from PHENIX [ 2 , 31] and the 
data in (b) and (c) are taken from STAR but have different pT ranges |32l. l33l| . 




Eigure 3: (Color online) Fitting results using the Tsallis distribution Eq. 
and p in p+p collisions at y/s = 900 GeV. Data are taken from (a) ALIGE 



for TT^, vr^, 

and (b) GMS 0. 


obtain 


(■E-^)hi<. -'4(1 + 


(13) 


where all the constants are absorbed into the new htting parameter A. n and T are the 
other two htting parameters. This form of Tsallis distribution was derived without resorting 
to thermodynamical description, which is different from Eq. (2). This also causes the power 
differences in Eqs. ([2]) and (ITSH . i.e. the power g/(l — q) = —nq in Eq. ([2]). Eq. flT3|) 
is equivalent to Eq. ([T]) but in a simpler form. We adopt it to do the systematic particle 
spectra analysis from p+p collisions. We will show that the approximation in Eq. (IT^ 


doesn’t change the results from Eq. 0 


been used by CMS Collaboration [7|, 


8l, 


later in this paper. We notice that Eq. (IT^ has 


28| and by Wong et al. in their recent paper 22]. 


The STAR Collaboration also applied a formula which is very close to Eq. (IT^ [29 1. 
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Figure 4: (Color online) Fitting results using the Tsallis distribution Eq. (I13p for vr^, K^, p and 
p in p+p collisions at y/s = 2760 GeV. Data are taken from CMS 0]. 

III. RESULTS 


We collect the spectra data for different particles with different pt ranges and different 
rapidity cuts from p+P collisions at y/s = 62.4, 200, 900, 2760 and 7000 GeV. A similar study 
has been done in ref. 19| assuming all the particles have the same temperatures while for the 
pions, which are produced the most in p+p collisions, the temperature is a free parameter. 
It is not necessary to require that the temperatures for all the particles are the same. If all 
the particles are produced at the same time and in thermal equilibrium, it is reasonable to 
have this constraint. If, for instance, the particles are produced at different times, i.e. in the 
framework of a cascade particle production mechanism, the particles will not have the same 
temperatures because of energy conservation. The particles produced at early time will have 
the higher temperature than the ones produced at later time. In this analysis, we allow the 
temperature to be a free parameter for different particles. This work is the natural extension 
of ref. 191]. Meanwhile, we include more data and a wider pt range. In our analysis, the 
pion case is exactly the same as in ref. 19| except the authors chose another form of the 
Tsallis distribution, see equation (7) in ref. 19||. For the other particles are different. Thus 
we cannot compare our results with ref. 19| directly. We will show this later in the paper. 

In Fig. [H we show our hts using the Tsallis distribution for different particles from the 
PHENIX Collaboration at y/s = 62.4 GeV. We can see that there is only little difference 
between the mesons and the corresponding anti-particles. The proton spectrum is over that 
of p because the colliding system is p+p and we have more protons. We can see that our 
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Figure 5: (Color online) Fitting results using the Tsallis distribution Eq. (11!^ for p 

and p in p+p collisions at ^/s = 7000 GeV. Data are taken from (a) ALICE [J] and (b) CMS p. 



Eigure 6: (Color online) Eitting results using the Tsallis distribution Eq. (1131) for charged particles 
in p+p collisions at ^/s = 900,2760 and 7000 GeV respectively. Data are taken from CMS [^,34]. 


fits with Eq. ffT3l) are good. The fitting parameters can be found in Table. [H 

Fig. [2] shows the hts for the identihed hadrons at y/s = 200 GeV with different pr ranges 
and rapidity cuts which are around mid-rapidity. As we can see, the particle spectra for 
kaons and protons start to split at high p^. Similar to Fig. [H our hts are excellent and the 
htting parameters can be found in Table. [B 

Similar to Figs. [H El we repeated the htting process for diherent particles produced at 
higher energies from LHC with diherent pt ranges and rapidity cuts. The htting results are 
showed in Figs. El 01 |5l As the collision energy gets higher, there is no diherence between 
the particle and its anti-particle as expected. Same to the htting results in Figs. [H El the 
htting quality in Figs. El 01 [5]are very good. We list the htting parameters in Table. [T] as 
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Figure 7: (Color online) Temperature T in the Tsallis distribution versus ^/s for (a) vr, (b)kaon 
and (c) proton in p+p collisions. 
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Figure 8: (Color online) The parameter n of the Tsallis distribution versus y/s for (a) tt, (b) kaon 
and (c) proton in p+p collisions. The lines are to guide the eyes. 


well. 


The prominent fitting power of the Tsallis distribntion is exhibited perfectly in Fig. | 6 l 
This is for charged particles. As we can see, the excellent htting can cover 15 orders of 
magnitnde np to 200 GeV/c for pt- This spectacular result was hrst showed by Wong et 
al 2l|. In Table. Uwe can see that the htting results are similar to the pions when we ht 
the spectra of the identihed hadrons separately at the same collision energy. Our htting 


211. This verihes that the 


results are consistent with the ones obtained by Wong et al. 
approximation we used in Eq. ffT^ is good. 

In Fig. [71 we show the parameter T of the Tsallis distribution obtained by htting spectra 
of particles at diherent ^/s in p+p collisions. Firstly, we can see that the parameter T 
depends on the htting pt range. Only the results for high px particles (black stars) are 
very much diherent from the ones including all the px particles at ^/s = 200 GeV. Secondly, 
T are almost the same for a particle and its anti-particle as it should be. Thirdly, T for 
pions produced at diherent collision energy i/s is almost a constant around 0.13 GeV. This 
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is consistent with the resnlts in a similar stndy for charged particles in p+p collisions 221 
and p+Pb collsions at ^/SNN = 5.02 TeV |35|. This may indicate that the pion production 
mechanism is similar in p+p and p+Pb collisions. T for kaons is a little bit smaller than the 
one for protons at the same collision energy y/s. But they have similar behaviors and increase 
with ^/s. We can interpret the results in the framework of a cascade particle production 
mechanism. Since the kaons and protons are heavier than pious, they are more likely to be 
produced at the beginning while pions can be produced at all times. Then their temperatures 
are higher than pions. When we include the results of the rj meson jd] and the multi-strange 
baryons and produced at y/s = 7000 GeV together, the trend is obvious. T is 

higher when the mass of the particle is heavier which is consistent with the cascade particle 
oroduction mechanism. Similar behaviors have been observed in p+Pb collisions as well 
35| . Finally, T for the identihed hadrons depends on the Tsallis distribution form used in 


the htting process. As we mentioned before, some authors argue that the identihed hadrons 


should have the same T 19|; others suggest otherwise m 


ise {2-5, 


15 


35|. Therefore we only 


can compare the results of pions from different versions of the Tsallis distribution. In this 
case, we also can put the results of charged particles from different versions of the Tsallis 
distribution in since T for charged particles is very similar to pions’s as showed in the 
previous discussion. We categorize Eqs. ([H[9l[13]) as type-A, Eqs. ([2l[3|) as type-B and Eqs. 
(ini [7|) as type-C Tsallis distribution to clarify the discussion. Type-A Tsallis distribution 


gives T ~ 0.13 GeV for pions and c 


. UE 


l arged p articles. This has been shown in Fig. 7] (left 

24. The 


21 


22|. Type-B gives T ~ 0.07 GeV 


14 


16 


panel) and can be found in refs 
extra term rriT in Eqs. (E] [3|) is responsible for this lower T, which is the main difference 


between Type-A anc 
gives T ~ 0.08 GeV 
gives T ~ 0.1 GeV. 

As discussed in refs. 


Type-B. Without it, Type-B will give similar T as Type-A. Type-C 


10| . Another type of Tsallis distribution, see equation (7) in ref. 19| . 


19 


20| . we can connect the parameter n in the power law of particle 


spectra at high px with the particle production process. We refer to [19|, l20| for details. Here 
we just briefly state the results. The parameter n = 2{na — 2) where Ua is the number of 
active participants in the point like scattering. Ua includes all the consituents of initial and 
hnal states in the scattering. If we assume that the basic scattering process at high px is 
qq —)■ qq, we have 2 constituents for initial and hnal states respectively, then = 2 + 2 = 4 
and n = 4. We can obtain higher values of n from the scattering at the hadron level, i.e. 
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p + meson —?■ p + meson, we have na = 5 + 5 = 10 and n = 16 which is the upper limit for n 
for pions and kaons. For proton, n can go higher up to 20 in the scattering process pp pp. 
This also sets the limits for the parameter q in the Tsallis distribution. Using Eq. ([S]), we 

Q G [1.05,1.25] for protons. The htting results of q 


obtain q G [1.063,1.25] for meso 

ns 

anc 


are between the limits in refs. 

10 , 

15, 

16] 


16| . In Fig. [HI we show the parameter n obtained 


by htting the spectra of particles with Tsallis distribution at different ^/s in p+p collisions. 
We can see that the values of n are well between those limits for all particles. This could 
give us some hints about the particle production mechanism in p+p collisions. We notice 
that the parameters n for pions and protons are decreasing with the collision energy which 
indicates that interactions are from hadron level to quark level. This picture is quite clear 
we can break the hadrons into quarks at higher collision energy. But if we look at the results 
from kaons, the n dereases hrst and then increases. Such a feature is not clear. 


IV. CONCLUSIONS 

We have made a thorough study of traverse momentum spectra of identihed particles 
produced in p+p collisions at RHIC and LHC energies with Tsallis distribution. A detailed 
analysis of the parameters T and n is also shown. T is not dependent on the beam energy 
for pions, while for kaons and protons it increases with increasing energy. Furthermore, 
we notice that T is higher for the particle whose mass is larger. This is probably due to 
the particle produced time. In the cascade particle production mechanism, this result is 
perfectly understandable. The behavior of n for kaons is not the same as pions and protons, 
which is related to the particle production process. From the properties of T and n, we get 
more information about particle production mechanism in p+p collision. We also wish more 
exciting results can be found in p+p collisions at 8 TeV. 
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Table I: The fitting parameters in Tsallis distribution Eq. (1131) for the particle spectra in p+p 
collisions. 


data source 

A/s(GeV) 

particle 

PT(GeV) 


T] 


A 

n 

r(GeV) 

PHENIX[30 

62.4 

TT^ 

0.614-6.734 

< 0.35 

0.128 

13.26 

0.135 

phenixH 

62.4 

Tr"*" 


0.35-2.85 


< 0.35 

0.104 

12.04 

0.132 




7T~ 


0.35-2.85 





0.111 

12.28 

0.130 




K+ 


0.45-1.95 





0.0035 

9.47 

0.160 




K- 


0.45-1.95 





0.0031 

11.38 

0.163 




P 


0.65-3.5 





0.0017 

11.95 

0.156 




V 


0.65-3.5' 





0.0017 

9.35 

0.126 

PHENIX[31 

200 

tt" 

0.616-18.93 

< 0.35 

0.218 

9.27 

0.118 

PHENIXI2J 

200 

Tr"*" 


0.35-2.95 


< 0.35 

0.174 

8.18 

0.114 




7T~ 


0.35-2.95 





0.148 

8.92 

0.123 




K+ 


0.45-1.95 





0.0065 

6.17 

0.136 




K- 


0.45-1.95 





0.0058 

7 

0.147 




P 


0.55-4.5 





0.0029 

8.32 

0.145 




P 


0.55-4.5' 





0.0022 

9.13 

0.152 

STARI32J 

200 

tt"*" 


0.35-9.0 


< 0.5 

5.219 

9.8 

0.129 




TT~ 


0.35-9.0 





5.077 

9.84 

0.130 




P 


0.468-6.5 

1 




0.073 

11.05 

0.181 




P 


0.468-6.5 





0.0629 

10.48 

0.175 

STARI33J 

200 


3.113-13.06 

< 0.5 

20.809 

8.87 

0.096 




TT~ 

'3.113-13.06' 




21.557 

9.04 

0.096 




K+ 

'3.113-13.06' 




7.666 

8.47 

0.083 




K- 

'3.113-13.06' 




6.123 

9.48 

0.092 




P 

'3.113-13.06' 




2.2782 

8.54 

0.0867 




P 

'3.113-13.06' 




2.2923 

9.21 

0.09 

ALICEna 

900 



0.11-2.5 


< 0.9 

5.333 

7.68 

0.125 




7T~ 


0.11-2.5 





5.279 

7.85 

0.126 




K+ 


0.225-2.3 





0.206 

5.81 

0.159 




K- 


0.225-2.3 





0.208 

6.17 

0.160 




P 


0.375-2.3 





0.0524 

7.05 

0.181 




P 


0.375-2.3 





0.0498 

8.08 

0.188 

ALICE 

4 


900 

TT^ 

0 

.495-5.818 

0 

0.26 

7.98 

0.135 

CMSI6J 

900 


0.125-1.175 

< 1.0 

6.139 

8.25 

0.134 




7r“ 

0.125-1.175' 




5.827 

8.18 

0.136 




K+ 

0.225-1.025' 




0.237 

6.54 

0.181 




K- 

0.225-1.025' 




0.235 

6.34 

0.179 




P 

0.375-1.675' 




0.0618 

7.02 

0.199 




P 

0.375-1.675' 




0.0598 

6.88 

0.195 

CMSI6J 

2760 

TT~^ 

0.125-1.175 

< 1.0 

7.702 

6.55 

0.128 




7T~ 

0.125-1.175' 




7.11 

7.21 

0.135 




K+ 

0.225-1.025' 




0.264 

10.84 

0.214 




K- 

0.225-1.025' 




0.258 

13.22 

0.218 




P 

0.375-1.675' 




0.0678 

6.39 

0.22 




P 

0.375-1.675' 




0.0632 

9.16 

0.238 

ALICE 

4 


7000 

tt" 

0.35-22.197 

0 

0.467 

6.84 

0.141 

CMSIC 

7000 

Tr"*" 

0.125-1.175 

< 1.0 

9.298 

6.4 

0.130 




7T~ 

0.125-1.175' 




8.815 

6.46 

0.133 




K+ 

0.225-1.025' 




0.302 

11.12 

0.228 




K- 

0.225-1.025' 




0.307 

7.57 

0.216 




P 

0.37^-1.675' 




0.0742 

5.93 

0.243 




P 

0.375-1.675' 




0.074 

5.97 

0.243 

CMS 9 


900 

charged 


0.5-31.2 


< 2.4 

15.938 

7.65 

0.128 
































































